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Abstract 



The low-energy excitations of filled Landau levels (LL's) of electrons in- 
volve promotion of a single electron from the topmost filled LL to the lowest 
empty LL. These are called excitons or collective modes. The incompress- 
ible fractional quantum Hall states are understood as filled LL's of composite 
fermions, and the low-energy neutral excitations are excitons of composite 
fermions. New techniques are developed to study large systems, which pro- 
vide detailed information about the dispersions of the composite fermion ex- 
citons. In particular, it is found that the interaction energy of the exciton is 

well described by the 'unprojected' composite fermion theory. 
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I. INTRODUCTION 



The quantum Hall effect (QHE) is one of the most fascinating discoveries in physics in 
recent years. It concerns electron transport in two dimensional systems in the presence of a 
high transverse magnetic field. It has been found that there are plateaus on which the Hall 
resistance is quantized fj]] at 



where / is either an integer or a simple fraction. For historical reasons, the observation of 

integer values of / is called the integer QHE (IQHE) and of fractional values of / is called the 

fractional QHE (FQHE). The longitudinal resistance in the plateau region is exponentially 

small, vanishing in the limit of temperature T — > 0. 

Not long after the discovery of the FQHE, Laughlin proposed the following wave 

function for the ground state at f — -: 

6 J (2m + 1) 

1>i/<Wri) = n(^-^) 2m+1 exph^EN 2 ] (2) 
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where Zj = Xj + iyj denotes the coordinates (xj , yj) of the jth electron, and the lengths are 

[he 

expressed in unit of the magnetic length Iq = \ — • This wave function was studied for small 

V e B 

systems and found to be quite accurate. Laughlin also wrote trial wave functions for the 
charged excitations (quasielectrons and quasiholes, collectively referred to as quasiparticles). 

Subsequently, Girvin, MacDonald and Platzman || (GMP) employed a single mode 
approximation (SMA) to investigate the neutral excitation of the Laughlin state, called the 
collective mode (CM) or exciton. The SMA takes the trial wave function of the collective 
excitation to be 

X S k MA = Vp^m+l) (3) 

where V is the lowest-Landau-level (LLL) projection operator and 

N 

P k = E ^ 3 (4) 



is the density wave operator with wave vector k, which is the Fourier transform of the density 

N 

p(r) = £ 6(r- r,) , (5) 

3=1 

for N electrons. The notation \ will be used exclusively for the wave function of the collective 
mode or exciton in this article; the superscript will denote which scheme is used, and the 
subscript will show the quantum number. GMP evaluated the energy of this wave function, 
given by 

< x SMA lH _ Eq1x SMA > 



< ^SMA^SMA > 



(6) 



as a function of k. In finite system studies, the SMA was found to work well in an intermedi- 
ate range of wave vectors near a minimum in the dispersion, termed the "roton" minimum, in 
analogy with Feynman's theory of superfluid 4 He |4j]. At wave vectors beyond the minimum, 
however, the SMA was not satisfactory. For large k, the energy of the true collective mode 
should be roughly independent of k, since it contains a far separated quasielectron-quasihole 
pair. The energy of xf MA , however, continues to grow with k. Laughlin attempted an alter- 
native description of the neutral collective mode in terms of a quasielectron-quasihole pair 
excitation g. This gave a reasonable description at large wave vectors, but did not produce 
the roton minimum. Thus, it was not possible to describe the full dispersion of the neutral 
excitation of / = l/(2m + 1) in a single theoretical framework. 

Much less was known about the collective modes of the other FQHE states. The SMA 
can be generalized to the other fractions as 

X S k MA = Pp^n/(2mn+l) , (7) 

TX 

where VWfamn+i) is the incompressible ground state at / = -. Since no good 

(2mn + 1) 

microscopic ansatz existed for ipn/(2mn+i), the SMA was tested in finite system studies, 
with the exact numerical ground state used for ip n /(2mn+i) ||- It, however, did not give a 
satisfactory description of the collective excitations JjJ. Finite system exact diagonalization 
studies also failed to provide a satisfactory overall picture for the collective excitations 



at general fractions. The results showed a significant amount of finite size effects, making 
an extrapolation to the thermodynamic limit difficult. 

In the last few years, there has been a resurgence of interest in this issue for two reasons. 
First, significant progress has been made on the experimental front. Pinczuk et al. || have 
measured the positions of the maxima and minima in the collective modes of the IQHE 
states, and recently also their dispersion in modulated density samples |§. Further, Raman 
scattering [|l(| and phonon absorption fill experiments have reported observation of the 



collective modes in the FQHE regime. 

Secondly, there now exists a new theoretical framework, called the composite fermion 
(CF) theory ||12|| , for describing all FQHE on an equal footing. The strongly correlated elec- 
trons are mapped on to weakly correlated composite fermions, which are electrons 'dressed' 
with 2m vortices of the many particle wave function. They may be loosely thought of as 
electrons carrying 2m flux quanta. The composite fermions experience a reduced effective 
magnetic field, and form energy levels, called CF-LL's or quasi-LL's, which are similar to the 
LL's of electrons in this reduced magnetic field. The FQHE of electrons is understood as the 
IQHE of composite fermions, and the FQHE ground state at / = n/ (2mn + 1) is interpreted 
as the state containing n filled CF-LL's. The Laughlin wave function, in particular, is in- 
terpreted as one filled CF-LL. This description of the incompressible FQHE states suggests 
that the low-energy neutral excitations can be obtained by promoting a single composite 
fermion to the next higher CF-LL (Fig. 1), i.e., they are excitons of composite fermions. 
Their wave functions can be constructed in analogy to the wave functions of the excitons of 
electrons for IQHE states, and exact-diagonalization studies for systems containing up to 9 
electrons have shown that they provide an extremely good quantitative description of the 
low-energy excitations of the FQHE states at 1/3, 2/5, and 3/7 for all values of k accessible 
in these finite size systems |T3|JT^1 . 

The excited composite fermion will be called a quasielectron of the FQHE ground state, 
and the hole left behind a quasihole [TJ|. The essential difference between the SMA and 
the CF wave functions is that the SMA creates an electron-hole pair, whereas the CF wave 



function creates a quasielectron-quasihole pair. We will see below that the SMA and the CF 

approaches are in general unrelated; the only exception is for the exciton at -, for 

(2m + 1) 

which the two become equivalent in the limit k — *■ 0. 

Despite the confirmation of the validity of the CF scheme in small system studies, it 

has not been possible to study large systems yet. The main difficulty is that the CF wave 

functions have a small amount of mixing with higher LL's to begin with, and must be 

projected on to the lowest LL to obtain good quantitative estimates in the large B limit. 

The projection operator is quite complicated and no general computational schemes exist 

for dealing with it for large systems. In the present study, we develop two techniques for 

computing the properties of CF excitons for large systems. The first one computes the energy 

of the LLL-projected wave function of the CF exciton of the state by generalizing a 

(2m + 1) 



method developed by Bonesteel [T^j. This method, however, is inapplicable to other FQHE 
states. In the second one, we show that the 'unprojected' wave functions themselves can 
be used to a great extent for the study of the exciton dispersion. We have computed the 
(self) energy of a single quasielectron or quasihole in the unprojected scheme, and found it 
to be off by a factor of two to three |T7| . However, the interaction between the quasielectron 
and quasihole participating in the exciton is found to be well described by the unprojected 
theory, except when they are very close. Therefore, the unprojected theory gives the exciton 
dispersion reasonably well up to an additive constant correcting for error in the self energies 
of the quasielectron and quasihole. We study large systems within the framework of the 
unprojected CF theory, and make qualitatively new and detailed predictions for various 
collective mode dispersions, which exhibit a rich structure with several minima. The extrema 
in the dispersion are of experimental relevance, since the CM density of states has peaks at 
the corresponding energies, which, as a result of a disorder-induced breakdown of the wave 
vector conservation, become observable in inelastic light scattering experiments ||. 

The analogy to the IQHE excitons is helpful in explaining the number of minima and 

Tl 

their positions. The relatively deep minima in the CM dispersion of the / 



(2mn + 1) 

FQHE state occur at approximately the same wave vector as those in the CM dispersion of 



the f = n IQHE state. The latter originate from an interplay between the density profiles of 
the quasielectron in the (n + l)th LL and the quasihole in the nth LL as they approach one 
another. Our study thus indicates that the wave function profile of the quasielectrons in the 

Tl 

f = FQHE state is similar to that in the f = n IQHE state, further deepening 

(2mn +1) 

the analogy between the FQHE and the IQHE implied by the CF theory. 

The dispersions for the collective modes have also been calculated in a Chern-Simons 
(CS) field theoretical scheme in several recent papers. Here, the composite fermions are 



modeled as electrons carrying (gauge) flux quanta |12],[l8l, which simulate vortices. The 
qualitative physical picture of composite fermions is obtained at a mean-field level, and a 
perturbation theory is carried out at the level of random-phase-approximation (RPA) |19| . 
Although the trial-wave-function and the CS schemes are both attempting to describe the 
same physics from different starting points, a precise relationship between them is not known 
at the moment. We will point out some similarities and differences between the results of 
the two methods in the specific context of the collective modes. 

The plan of the paper is as follows. In Section II, we write the wave function of the CF 
exciton, and show that, for / = l/(2m+ 1), it is equal to SMA wave function in the limit of 
small wave vectors. Section III gives the results of our study of the projected CF excitons, 
both in small systems, where a comparison with the exact results is possible, and in large 
systems, where Monte Carlo techniques allow the computation of its dispersion for 



the 1/3 state. In Section IV, we study the interaction energy of the unprojected CF-exciton 
wave functions. It is shown that it provides a reasonably good quantitative description of the 
interaction between the quasielectron and the quasihole. Section V considers the analogy 
between the excitons of composite fermions, relevant for the FQHE, and those of electrons, 
relevant for the IQHE. This provides insight into the results of the previous sections, and 
also sheds light on the origin of the minima in the FQHE exciton dispersions. Section VI 
investigates the kinetic energies of the unprojected wave functions. It is shown to be less 
than the cyclotron energy for finite k, approaching the cyclotron energy as k — > 0. This 
is reminiscent of the Kohn's theorem, a weaker form of which is given for a class of wave 



functions which are not eigenstates of any known Hamiltonian. The paper is concluded in 
Section VII. A short report of parts of this paper has been published elesewhere |Z2| . 

For future reference, here is a list of notations we will be using below for various wave 
functions and energies. 

$ n : Slater determinant wave function of n filled LL's of electrons. 

(§n P CF )'- Projected (unprojected) wave function of n filled LL's of composite fermions; 
this corresponds to the ground state of electrons at / = nj (2mn + 1). The superscript 'UP' 
denotes 'unprojected'. 

X {x UP ~ CF ) : Projected (unprojected) wave function of the CF exciton. 

AV ex : Exact energy of the exciton (measured relative to the exact ground state) from exact 

diagonalization in the lowest Landau level. 

AV p : Interaction energy of x relative to <& CF . 

AV (AK): Interaction (kinetic) energy of X UF ~ CF relative to §> UP ~ CF . 



II. COMPOSITE FERMION EXCITONS 



The basic principle of the CF theory |12| is that, in a range of filling factors, the electrons 



in the lowest LL find it energetically favorable to capture an even number of vortices of the 
many particle wave function. The bound state of an electron and the vortices behaves as a 
particle, called composite fermion. The vortices produce phases as the composite fermions 
move around, which partly cancel the Aharonov-Bohm phases originating from the external 
magnetic field, and, as a result, the composite fermions experience an effective magnetic 
field given by 

B* = B- 2mp0 o , (8) 

where B is the external field, 2m is the number of vortices carried by composite fermions, 
4>o = hc/e is the flux quantum, and p is the electron (or CF) density. The residual interaction 
between the composite fermions is weak, and the strongly correlated liquid of electrons maps 
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into a weakly interacting gas of composite fermions. An effective single-particle description 
of the electron state thus becomes possible in terms of composite fermions. The energy 
levels of composite fermions are analogous to the LL's of non-interacting electrons in this 
weaker magnetic field, and are called quasi- or CF-LL's. Defining the CF filling factor as 
v* = p<po/B*, in analogy to the electron filling factor v = p(f) /B, the above equation can 
also be expressed as 

V ~ (2mu* + 1) ' ^ ' 

The IQHE of composite fermions at v* = n manifests as the FQHE of electrons at 

n 



v 



{2mn + 1) 

The trial wave functions for the CF states are constructed as follows. For simplicity, 

we confine our discussion to the special filling factors v = -, which corresponds to 

(2mn + 1) 

the CF filling factor v* = n. Let us denote the ground state of non-interacting electrons at 
v* = n by $„. The wave function of n filled CF-LL's is obtained by attaching 2m vortices 
to each electron in the state $ n , which amounts to a multiplication by the Jastrow factor 

D m = n( Zj - z k f m , (io) 

j<k 

which produces 

< P - CF =m^-^f m ^n. (11) 

j<k 

This wave function is not strictly in the lowest LL (although it is largely so ||23||), and must 
be projected on to the LLL to obtain a wave function appropriate to the limit B — > oo: 

*™ = V^ n p - CF . (12) 

Ti 

describes the electron ground state at v — -. It has been tested for several 

(2mn + 1) 



FQHE states in small system exact calculations and found to be extremely accurate fll"3| , |T4 

Armed with this analogy between the IQHE and the FQHE states, we are now in a 
position to write a wave function for the CF exciton (or the CF collective mode). It is 
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obtained from the exciton of the v — n IQHE state precisely in the same manner as the 
FQHE ground state is obtained from the v = n IQHE ground state, i.e., by multiplication 
by the Jastrow factor, followed by LLL projection: 

X^ = -P X ^ P - CF , (13) 

where 

X UP-CF = n(z . _ Zk) 2m ^n+r^ (M) 
j<k 

Here, p™^ n+1 Q n is the exciton of the v = n state. In this wave function, a single electron has 

been excited from the vP^ 1 LL of $ n to the (n + l)^ 1 LL, creating an exciton whose kinetic 

energy is equal to the cyclotron energy at the reduced magnetic field, Tiui* = — - p3 . 

(2mn + 1) — 

By analogy, x CF contains a single excited composite fermion in the {n + l)" 1 CF-LL and a 

th 

hole left behind in the n CF-LL, as shown schematically in Fig. 1. 

The SMA wave function for the collective mode at v — r has a CF-type inter- 

(2m +1) JF 

pretation, as it is obtained from after multiplication by the Jastrow factor and then 
projecting the product on to the lowest LL. We now show that the CF and the SMA wave 
functions become identical in the limit of k — > for this mode. This follows because, in 
general, in the limit k — > 0, 

P k ^n = C; +1 $n ■ (15) 

To see this, we express p k and p™^ n+1 in the second quantized form. Denote the single 
particle states by \p,s >, where p is the wave vector quantum number and s = 0, 1, ... is the 
LL index (note that the topmost occupied LL of $ n is s = n — 1). The explicit form of the 
wave function in the Landau gauge is given by 

P , S = (2vr2 s S !v^r 1/2 exp[ipy -\(x + p) 2 ] H s {x + p) , (16) 

where H s is the Hermite polynomial, and the magnetic length has been set equal to unity. 
In the second quantized notation, the density operator is given by 
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P k = E al, s ,a ps <p's'\e iky \ P s> (17) 

p,p',s,s' 

where the y-axis has been chosen parallel to the wave vector k. In the present case, we have 
s < n and s' > n since p k is applied to a state with n filled LL's. Aside from an overall 
normalization factor, it is given by, for s' > s, 

Pk = ^ s+s 's\s'\)- ll2 e- k2 lH- s s\k s '- s Li- s (e/2) J dpal +k>s ,a ps . (18) 

s,s' 

The operator p™^ n+1 is given by the term with s = n — 1 and s' = n. Eq. fll5|) follows in 
the limit k — > 0. For fractions other than v = l/(2m + 1), -^ SMA does not have a CF-type 
interpretation, and y° F and yf MA are different at all k. 



III. PROJECTED WAVE FUNCTION 



This section investigates the projected CF wave function. The spherical geometry P5|J21] 



is used in all our calculations, where N electrons move on the surface of a sphere under the 
influence of a radial magnetic field. The flux through the surface is 2q<po, where 2q is 
an integer. The single particle eigenstates are called monopole harmonics, Y qj ^ z , where 
^ = I? I? \q\ + 1) •••is the single particle angular momentum, and £ z = ±£, ±£ — 1, ... is the 
z-component of i. The total orbital angular momentum L of the many body system is related 
to the wave vector of the planar geometry by kl = — . In the spherical geometry the 

energy separation between the rfi 1 and the (n + 1)^ LL is given by ( 1 H — ] — - — , which 

\ qj mc 

reduces to the usual cyclotron energy in the limit q — > oo. We find it more convenient to use 

as our unit of the kinetic energy hu c = ( 1 H — ] , the separation between the lowest 

\ q J mc 

two LL's. 

If the highest occupied LL shell in $ n has angular momentum £, then the IQHE exciton 

t h 

has a single excited electron in the (n + 1) LL, with angular momentum t + 1, and the 
th 

hole in the n LL, with angular momentum £. The Slater determinant basis states can be 
denoted by 

\et,e h z > 
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where t e z (i 1 ^) is the z-component of the angular momentum of the electron in the (n + 1) 

f h 

LL (hole in the n LL). The allowed values of L for the exciton are L = 1, 2 ...21 + 1, with 
precisely one multiplet at each L, containing 2L + 1 degenerate states. The k = in the 
planar geometry corresponds to L = 1 in the spherical geometry, which is the smallest L 
possible for the exciton. With no loss of generality, we restrict our calculations below to the 
sector with L z = 0, with the understanding that each state in this sector represents 2L + 1 
degenerate states of the full Hilbert space. The exciton wave function p™~^ n+l Q n is a unique 

linear superposition of the Slater determinant basis states, given by (with L z = 0) 

i 

P r n+1< ^= <e + i,e g -,e,-e s \L,a> |4,-4>, (19) 

The wave function of the CF exciton at the corresponding L, x? F \ is obtained by multipli- 
cation by the Jastrow factor, as in Eq. ([13]). The Jastrow factor in the spherical geometry is 
$1, where <3>i is the wave function of the lowest filled LL, multiplication by which does not 
change the L of the state. Clearly, x° h F also does not contain any adjustable parameters. 

We have computed the exciton energy using the projected CF wave function x CF \ denoted 
by AV P , for finite systems. It is measured relative to the projected CF ground state energy. 
Our brute force projection method (see Ref. JHJ for details) allows us to carry out the 
projection for general states only for up to 8-9 electrons, which severely limits the size of 
the systems for which AV P may be evaluated. 

Figs. 2 and 3 give the actual Coulomb energies, AV ex , obtained by an exact numerical 
diagonalization of the Hamiltonian, along with the energies of the projected CF wave func- 
tions. First note that the range over which the exciton extends in the exact diagonalization 
studies agrees with that predicted by the CF theory. (A state exists at L = 1 for the IQHE 
exciton, but it is annihilated by the projection operator while constructing the wave function 
for the CF exciton.) AV P is within a few percent of the exact energy AV ex . 

For the exciton at v — 1/3, we have obtained AV P for large systems using variational 
Monte Carlo, which relies on the special feature of x UP ~ CF that it contains no more than one 
electron in the second LL, and none at all in the higher LL's. This follows since the IQHE 
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exciton wave function, p 1 ~ r2 <&\, has only one electron in the second LL, and multiplication by 
the Jastrow factor does not promote electrons to higher LL's. The projected wave function 
can then be written as [[U] 



X CF oc (K - hco c ) X UP - CF (20) 
1 eB 

where, we remind the reader that huo c = (1 H — )h — . The kinetic energy operator K, in 

q mc 

standard notation is given by 

1 ( d ( . „ d \\ I 



sin ^^r)) + -^(gcos^-mj) 2 



N 

- -hu c (21) 



sin^- \d9j \ J d6,j) ) sin 2 ^- 

where our choice of the gauge (or, 'section' f25fl ) is defined by id<f>Y q ^ z = (q — £ z )Y q ^e z , 
and the kinetic energy is measured relative to the LLL. Each Slater determinant is updated 
at every step of the Monte Carlo following the method of Ref. [^TJ, which takes 0(N 2 ) 
operations. The computation of the energy is rather involved and time consuming not only 
due to the presence of derivatives in the expression of K , but also from the fact that x UP ~ CF 
is a linear superposition of ~ O(N) Slater determinants, necessitating 0(N 3 ) operations for 
updating the determinants at each step. It should also be mentioned that the computation 
of the interaction energy of x° F is much more stable when \x UP ~ CF \ 2 1 rather than |x CF | 2 , 
is chosen as the Monte Carlo weight function. 

The CF-exciton energy is shown as a function of L for a 20 electron system in Fig. 4. 
Computation of each point takes approximately 5-6 days of computer time on the IBM 
RISC 6000 workstation. The k — > (L = 1) energy, 0.15e 2 /eZo, is consistent with the SMA 
prediction. There is a deep minimum at the expected position of kl w 1.4. From Fig. 5, 
the estimate for the thermodynamic value of the energy at the minimum is 0.063(3)e 2 /e/o, 
which should be compared to the SMA prediction 0.078e 2 /e/ [[§• We note here that the 
same wave function for the ground state (the Laughlin wave function) is used in both the 
CF and the SMA calculations. 

Two additional minima are clearly visible at Hq ~ 2.7 and Hq ~ 3.5. Are they real? 
We believe so. A direct confirmation of the genuineness of the minimum at Hq « 2.7 is 
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seen in the nine-electron exact diagonalization calculation of Fano et al. |27|| , reproduced in 
Fig. 6(c). (The exact diagonalization systems are too small to see the third minimum.) We 
note that there is no principle that rules out the existence of more than one minimum in the 
dispersion of the 1/3 exciton; the structure here arises simply from an interplay between the 
structures in the density profiles of thequasielectron and quasihole, as the distance between 
them is varied. This will be discussed later in greater detail. 

We note that the exciton energy predicted by the projected CF theory may not be very 
accurate in the limit k — > 0, as indicated by the following considerations, (i) The LLL 
projection of the unprojected state becomes small and vanishes as k —>■ (i.e., at L = 1), 
suggesting that the act of projection may not be as inoccuous here as at larger k. (ii) It was 
found in small system calculations that x CF does not exhaust the oscillator strength at small 



wave vectors to the same extent as it does at larger k [|13||14| . This was attributed to the fact 
that the exciton branch comes quite close to the higher bands at small k; the exact state 
mixes with higher band states to reduce its energy, whereas any such mixing is neglected in 
the CF approach. (For 2/5 and 3/7, the CF description continues to work up to the smallest 
wave vectors available in the finite size studies.) (iii) In the SMA theory, the k — > energy 
of the 1/3 collective mode was close to twice the energy at the minimum. Our somewhat 
improved estimate shows that the ratio between the two energies is definitely greater than 
two, approximately 2.4. This suggests that the actual lowest energy excitation at small k 
may contain two pairs of CF-excitons — these will have a more complicated quadrupolar 
structure [j3p8|. A better understanding of the nature of the low-energy excitations in the 
small-fc limit would require further investigation. 

For other fractions, <fr Up - CF ane x UP ~ CF have a finite fraction of electrons in higher 
LL's, and the above method for the LLL projection is not useful. Below, we use a different 
approach which, though less accurate, can be applied to the excitons of other FQHE states 
as well. 
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IV. UNPROJECTED CF EXCITON: INTERACTION ENERGY 



It has been known that the unprotected CF wave functions already reside mostly in the 



lowest LL [23], with very small kinetic energies per particle. It is natural to ask if they 
can be used, without projection, for the computation of various quantities. This has been 
addressed for the charge gap (relevant for transport), equal to the energy required to create 
a far separated quasielectron-quasihole pair, which is nothing but an exciton at a large k. It 
was found that the gap computed with the unprojected wave functions was smaller than the 
projected gap by as much as a factor of 2-3 JI7| . This may suggest that the unprojected 
wave functions are not relevant for a quantitative understanding. 

Fortunately, this turns out not to be the case. We will see in this section that the 
interaction energy of x UP is close to the energy of x CF \ up to an overall additive constant, 
except at very small k. In other words, if the interaction energy of x UP is called AV 
(measured relative to the unprojected CF ground state), then AV^ + T provides a good 
approximation to the exciton dispersion, except at small k. The constant V can be fixed by 
requiring the large k limit to equal the transport gap (which is experimentally measurable). 
The underlying physics will be discussed in the next section. 

We first test this assertion in finite system calculations. The energy AV^ + V (with a 
suitable choice of V) is shown in Figs. 6 and 7 for finite systems, where AV^ is computed by 
variational Monte Carlo. It indeed captures the essential features of AV ex ; in particular, it 
obtains correctly the minima and maxima. It also provides a good quantitative approxima- 
tion for the exciton energy: it is quite accurate for 1/3, and reasonably accurate (to within 
15-20%) for 2/5 and 3/7 (for which the agreement improves with the system size). 

The advantages of working with the unprojected CF wave functions are that a treatment 
of large systems becomes possible, and that all FQHE states can be considered. Fig. 4 gives 
the dispersion (AV + T, where T = 0.064e 2 /e/o is obtained from the extrapolation of finite 
system values) for a 20-electron system for the 1/3 FQHE state. A lack of any significant 
size dependence for systems with slightly larger N shows that these results are close to 



14 



the thermodynamic limit. Let us first concentrate on the range klo > 0.5. Here, AV + T 
provides a good approximation for AV p . A deep minimum appears at the expected position 
klo ~ 1-4. The other two minima at klo ~ 2.7 and klo ~ 3.5 are also seen in AV, although 
they are less pronounced than in AV P . AV also underestimates the size of the structure in 
Fig. 7, which seems to be a general feature. 

AV + T is shown in Fig. 8 for 1/3, 2/5 and 3/7 for large systems, with V suitably chosen 
to produce a reasonable large- k limit. The positions of the two deep minima for the 2/5 
exciton agree well with those found in the exact diagonalization results of |J, as also does 
the feature that the second minimum is deeper. The relatively complicated structure in 
the dispersion clarifies why the small system calculations are unable to provide a coherent 
picture. 

A curious feature of the exciton dispersion in Fig. 4 is that AV bends downward at small 
wave vectors (klo < 0.5 for 1/3). The unprojected scheme is not trustworthy here, since LL 
mixing becomes a more serious problem here. (Remember, at L — 1, x UP ~ CF has a zero 
projection on the lowest LL ||T^[). Indeed, AV P for 1/3 shows no such bending. 



V. ANALOGY TO IQHE 

An insight into the above results can be gained from a comparison with the exciton of 



the v* = n IQHE state [24 



A. positions of minima 

First we consider the number of minima and their positions. Fig. 9 shows the exciton 
dispersion of the IQHE states at v* — 1, 2 and 3. These have one, two and three minima 
(some of these are relative to a rising background, and appear as inflection points), respec- 
tively, which corresponds exactly to the number of strong minima in the exciton dispersions 
for 1/3, 2/5 and 3/7. To pursue this analogy further, we have shown by arrows the positions 
(L) of the minima/inflection points of the corresponding IQHE state in the dispersions of 
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Figs. 2- 7 (in each case, the IQHE state with the same N has been used to determine the 
minima). For the 1/3 state, this obtains the L of the minimum off by one unit (which is 
an insignificant error in the thermodynamic limit) while for 2/5 and 3/7, it predicts the 
positions of minima correctly. Given that L = kl 0y /q, and Iq 1 oc \fB oc y/q, this implies 
that the minima occur at the same k for the FQHE and IQHE states related by the CF 
theory. 

B. Interaction energy of CF exciton 

The exciton contains a quasihole and a quasielectron. For sufficiently large k, when the 
quasiparticles are far apart, the energy of the CF exciton can be thought of as having three 
contributions: the self energy of the quasielectron (which is the energy to create an isolated 
quasielectron), the self energy of the quasihole, and the interaction energy between them. 
The interaction energy is a monotonic function of the distance between the quasiparticles, 
proportional to 1/r, when they are far apart. As they approach one another, there is a 
non- monotonic, oscillatory contribution to the interaction energy originating from the fact 
that the the quasiparticles are not point-like objects but have a finite extent, with oscillatory 
structure in their density profiles. At very small k, when the size of the exciton becomes of 
the order of the size of a single quasiparticle, the distinction between the self and interaction 
energies becomes nebulous, and the above intuitive picture breaks down. 

As mentioned earlier, the unprojected CF theory is not able to predict the self energies 
of the quasiparticles quantitatively. However, these do not have any k dependence, and 
can be corrected by an additive constant Y. The results of the previous section show that 
the unprojected CF approach does a good job of estimating the interaction energy of the 
quasihole and quasielectron, with is responsible for the structure in the exciton dispersion. 

So long as the interaction energy is small compared to the self energy, the densities of 
the constituent quasiparticles are not much affected. The structure in the interaction energy 
then arises largely from the structure in the density profiles of the individual quasiparticles. 
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This motivates us to ask if the unprojected CF theory gives a reasonable description 

of the density profiles of the quasiparticles. To this end, we first investigate if the FQHE 

quasiparticles have any relation to the IQHE quasiparticles. The simplest case is of the 

1 

quasielectron and quasihole of the state. These are related, respectively, to the 

(2m + 1) 

quasielectron and quasihole of the v* — 1 state, i.e., to states with the LLL completely 
occupied except for a hole, and the LLL completely occupied plus an additional electron in 
the second LL. The associated densities are given by 



e"^ 2 , (22) 
Po 



and 

p(r) 



1 + , (23) 



Po 2 

respectively, where po = (2vr/g) _1 is the density of the filled LL. These are plotted in Fig. 10. 
The density profiles are different simply because the hole is in the LLL, while the additional 
electron resides in the second LL. The analogy to the IQHE thus provides a natural expla- 
nation of the asymmetry between the quasihole and quasielectron. The charge densities of 
the quasihole and the quasielectron of the Laughlin 1/3 state |JU| are similar to those in 



Fig. 10. In particular, the quasielectron of the 1/3 state has a smoke- ring shape with a dip 
at the origin. The maximum of the quasielectron of the v* — 1 state is at r = \/2 l . This 
predicts a maximum for the 1/3 quasielectron at r = V6 lo (since the lo's at v — 1 and 
v — 1/3 differ by a factor of v^) , which is in excellent agreement with r « lo found 
in earlier studies [|30|. The value of p/po at the origin is 1.0 and that at the maximum is 
1 + e _1 = 1.37, which should be compared to p/po ~ 0.93 and p/po ~ 1.26 (where p now 
is the electron density of the 1/3 state), respectively, for the 1/3 quasielectron (taken from 
Morf and Halperin [|30|). The quasiholes at v = 1/3 and v = 1 also have rather similar 
shapes, except near the origin. The similarity between the quasiparticles of a FQHE state, 
which are completely determined by interactions between electrons confined to the lowest 
LL, with the quasiparticles of the corresponding IQHE state, governed by the LL physics, 
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is quite impressive. This appearance of the higher-LL-like structure within the lowest LL is 
an explicit illustration of the basic principles of the CF theory. The wave functions of the 
quasiparticles of other FQHE states are expected to have more complicated profiles, since 
the radial part of the single particle wave function in higher LL's has several nodes. 

Multiplication by the Jastrow factor does not change the main features of the IQHE 
quasiparticles, although it introduces additional weaker oscillations. It was seen in Ref. [16| 



that x UP CF gives a good description of the density profile even in the region of overlap 
between the quasielectron and quasihole. 

C. Small k limit 

Another feature of the unprojected results is that AV^ decreases for small k (apparently 
vanishing as k — > 0). This can again be understood by analogy to the IQHE collective 
modes, for which also the interaction energy vanishes in the limit k — > 0, as guaranteed by 



the Kohn's theorem \Ttf[. 



The analogy between the FQHE and the IQHE excitons is, of course, not true in every 
detail. In particular, it applies only to the strong structure. There are weaker minima in the 
CF exciton dispersion, which have no analog in the dispersion of the corresponding IQHE 
exciton. These appear only after multiplication by the Jastrow factor (with or without the 
LLL projection), and thus refer to very fine features beyond the mean-field theory. These 
are related to weaker oscillations in the density profiles of the FQHE-quasiparticles that 
are not seen for the corresponding IQHE-quasiparticles. Also, the relative strengths of the 
strong minima do not necessarily correspond to those in the IQHE dispersion. For example, 
at v* = 2, the second minimum is much weaker than the first one (Fig. 9), unlike at v = 2/5 
where the first one is weaker (Fig. 8). 

The principal minimum for the 1/3 collective mode has been interpreted as a precursor 
of the Wigner crystal instability ||. The above interpretation is cast in terms of the IQHE 
of composite fermions. A reconciliation between the two should prove interesting. 
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VI. UNPROJECTED CF MODE: KINETIC ENERGY 



Further insight into some features of the above results can be gained by considering the 
kinetic energy of the various unprojected wave functions. We emphasize that it has no direct 
relevance to the actual energy of the collective mode (at least in the limit of large B), except 
that it must be sufficiently small in order for the LLL projection to retain the important 
correlations of the CF theory. 

The kinetic energy of the IQHE exciton (measured relative to the ground state) is pre- 
cisely hu*. We consider here the kinetic energy AK of the CF exciton, defined as 

AK <X UP - CF \K-E \ X ^> 

< yUP-CFVyVP-CF > ' 

where 

< $up-cf\ k \$up-cf > 

E = ~ InP-CFi^TIP-CF Z ' ( 25 ) 



and K is the kinetic energy operator. The kinetic energy per electron for the ground state, 



Eq/N, is known to be small [23 



Before coming to numerical calculations, we mention an exact result. It was proved 
earlier that the rigorous upper limit for AK for v = I /(2m + 1) is the cyclotron energy, 
since x UP cannot have more than one electron in the second LL. It has also been shown 
earlier that at L = 1, the wave function has zero projection on to the lowest LL, indicating 
that it has at least one electron in the second LL f/||l3]. Taken together, one concludes that 
X UP ~ CF has precisely one electron in the second LL at L = 1, i.e., its energy approaches the 
cyclotron energy in the limit k — > 0. For other fractions, no exact results are known. 

AK cannot be evaluated analytically in general. We have computed it by two methods. 
For small systems, we expand x UP ~ CF an d § UP ~ CF in the Slater determinant basis. For 
X UP ~ CF i there are a total of 18,814 basis states for N = 7 at v = 1/3 (in the L z = 
subspace, allowing a maximum of one electron in the second LL), and 4,448,884 basis states 
for iV = 8 at v = 2/5 (in the L z = subspace, allowing for a maximum of four electrons in 
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the second and one electron in the third LL). The amplitudes of various basis states can be 
obtained by an elaborate book-keeping. Despite the large basis sizes, we are able to evaluate 
AK exactly (numerically) for these systems (Fig. 11), since it is not necessary to store the 
amplitudes of all states. In these calculations, we put by hand the LL spacing to be the 
same for all consecutive LL's. As expected, AK of the 1/3 state at L = 1 is exactly equal to 
the cyclotron energy. For 2/5 also, AK at L = 1 is close to hu c as seen in Fig. 11(b). The 
actual values for AK/hu c at L = 1 for 4, 6, and 8 electrons are 1.0219, 1.0091, and 1.0050, 
respectively, approaching unity approximately as 1 + 0.9iV -5 / 2 . 

We have studied bigger systems by Monte Carlo, and the results are shown in Fig. 12. 
AK is small, less than the cyclotron energy in the entire k range for all FQHE states. It 
approaches Tiuj c as k — > 0, in accordance with the Kohn's theorem (below). At large k, its 
value is equal to the kinetic energy of a far separated quasielectron-quasihole pair. For the 
1/3 state, the large-A; limit is equal to 0.16 hu> c , which is also the kinetic energy of an isolated 
quasielectron |23] (since the quasihole is strictly in the lowest LL in this case). For the 2/5 



and 3/7 states, AK shows 2 and 3 shallow minima respectively which occur at approximately 
the same positions as the minima in AV UP ~ CF , demonstrating that AK also shares some 
qualitative features of the true exciton dispersion. Also note that AK increases at small 
k, which sheds light on why the results obtained in previous sections become questionable 
here. 

As discussed earlier, the SMA wave function for the collective mode at 1/3 is obtained 
by CF transformation from the state p k $i- The kinetic energy of this wave function is given 
by 

<$ 1 \p{(K-E ) Pk \$ 1 > 



< $i|plp fc |$i > 
Substituting from Eq. (|I8"D, this is equal to 

Em=i mS(m 



(26) 



a*s(m) (27) 

with S(m) = 2 m m\k 2m , where we have used the fact that L™(x) = 1 for arbitrary x. This 
yields the kinetic energy of p fe $i to be 
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It approaches the cyclotron energy in the limit k — » 0, as expected from Kohn's theorem. As 
k increases, however, the kinetic energy of increases, as shown in Fig. 13 (solid line), 
indicating an increasing occupation of higher LL's. Clearly, it is not an appropriate wave 
function for the exciton of the v = 1 IQHE state. The SMA wave function thus has a large 
occupation of higher CF LL's at large k. Fig. 13 also shows (squares) the kinetic energy of 
the 'unprojected SMA wave function', Y[j<k( z j ~ z k) 2 p k &i, measured relative to the kinetic 
energy of the state Ylj^izj — Zk) 2 $i, for 20 electrons. It approaches Tiuj c as k —>■ for the 
same reason as AK (also see, He et al. 0), but increases rapidly beyond kl ~ 1.5. 

Before closing the Section, we state a generalization of the Kohn's theorem, which sheds 
light on the feature seen above that the kinetic energy of the unprojected wave functions 
approaches huj c in the limit of k —>• for all FQHE states. The Kohn's theorem states that 
for any given eigenstate with eigenvalue E, is also an eigenstate with eigenvalue 
E + hu c , with 



N 



oJ = XX , (29) 

t=i 

where a\ is the LL raising operator. The unprojected CF wave functions, however, are not 
eigenstates of the Hamiltonian (in particular of the kinetic energy). A weaker form of Kohn's 
theorem applies to a class of wave functions which are not eigenstates of the Hamiltonian. 
More specifically, 

° ' tT , t ' ° = fcwc, 30 

where H = K + V and E =< ^f\H\^f > / < >, provided that \I/ can be written in a 
factorized form 

* = $ Ccom (31) 

where (i) $ does not depend on the center-of-mass (COM) coordinates, (ii) (com depends 
only on the COM coordinates, and (iii) (com is an eigenstate of the COM kinetic energy. 

21 



All these properties are satisfied by any eigenstate. This generalized Kohn's theorem can 
be proven straightforwardly by transforming to the center of mass and relative coordinates, 
and noting that there is no mixing between the COM and relative coordinates in the kinetic 
energy K, which is a sum of two parts, one containing only the COM coordinate and the 
other containing only the relative coordinates. The operators a and a® are recognized as the 
LL raising and lowering operators for the COM coordinate. The theorem then follows since 
the application of aj changes the kinetic energy by precisely hu c , but leaves the interaction 
energy unchanged (since it is independent of the COM coordinate). 

The unprojected CF wave function for the incompressible FQHE state at v — n/(2n + 1) 
is given by 

1 N 

D m K exp[--5>| 2 ] , (32) 
4 i=i 

where the prime on $^ denotes that the exponential has been factored out. To see that 
it has the same form as the wave functions in Eq. flSOj), note that the exponential can be 
expressed as 

1 X- 1 N-l N 

ex P[-7(H(^ -^o)« -^o) + I ( z i ~ Vo)\ 2 )} exp[-— W q] , (33) 
4 i=i i=i 4 

where r/o = A^ -1 J2f=i z j is the COM coordinate. The factor D m &' n is independent of the 

COM coordinates, since it is invariant under the replacement Zj — > Zj — r/ ). The COM 

coordinate thus factors out; in fact, (com — exp[— ^-^o^o] is a ground state of the COM 

kinetic energy. 

These considerations provide a systematic way of constructing new trial states at the 
same energy and also at energies differing by Hu c , by using the angular momentum and LL 
raising operators of the COM coordinate. This is related to the fact, seen routinely in finite 
size studies, that once an eigenenergy appears at any angular momentum value, it appears 
at all larger angular momenta as well |32|.|33| . 
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VII. CONCLUSION 



In short, the following unified picture becomes possible in the CF framework: the incom- 
pressible ground states contain an integer number of filled CF-LL's and the excited states 
are obtained by exciting composite fermions to higher CF-LL's. Not surprisingly, the same 
physics provides the key to understanding both the ground state and its excitations. 

We have carried out a detailed study of the neutral excitations of the 1/3, 2/5 and 3/7 
states and developed techniques that allow us to deal with fairly large systems. Predictions 
have been made for the rich structure of the collective modes of various FQHE states. The 

Tl 

dispersion of the r state has n strong minima, and several additional weaker minima. 

(2n + 1) 

The strong minima can be explained by analogy to the IQHE, and appear at the same wave 
vector as the minima in the collective modes of the IQHE states. The weaker minima, on 
the other hand, have no analog in the FQHE. The existence of these structures has been 
corroborated, whenever possible, with the help of exact diagonalization study. 

There is qualitative agreement between the results of our calculations and those of the 
Chern-Simons theory. The latter also obtains flat dispersion at large k, several minima, and 
finite energy in the limit k —>■ 0. A quantitative comparison is not as good, however. A 
crucial feature of the CS scheme is that, at small k, the mode derived from the n —>• n + 1 
IQHE mode is pushed up to the cyclotron energy due to the RPA screening. This may be 
related to the feature discovered in Section VI that AK — > Tiuj c as k — > for the unprojected 
CF wave function. A better understanding of these and other issues will require further 
work. 

Several effects left out in the above study must be incorporated before a comparison with 
experiment may be made. Modification in the Coulomb interaction due to the finite width 
of the quantum well [34|, LL mixing [35], and disorder |]5|j36[1 are all known to change the 
numerical values of the excitation energies. We believe that a good first approximation for 
the experimental CM dispersion can be obtained by AV + V, with a choice of F that makes 
the large-A; limit equal to the experimental transport gap. 
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FIG. 1. Schematic view of a CF exciton. Each filled circle with two arrows depicts a composite 
fermion as an electron carrying two flux quanta. All composite fermions on a given line are in 
the same CF Landau level. A CF exciton of the n filled CF-LL state [which corresponds to 
v = n/(2n + I) FQHE state of electrons] is obtained by promoting a single composite fermion from 
the topmost occupied [nth] CF Landau level to the lowest unoccupied [(n+ l)th] CF Landau level. 
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FIG. 2. The exact energies AV ex (dashes) and projected CF energies AV P (triangles) are shown 
for the exciton of the 1/3 state for systems of 6 and 7 particles. All interaction energies here and in 
the following curves are in units of e 2 /elo, where Iq is the magnetic length, and e is the background 
dielectric constant. The arrows here and in the following figures indicate the angular momenta at 
which minima or inflection points occur in the exciton dispersion for the corresponding filled LL 
IQHE state (y = 1 state in this figure) for the same N. The energy of the ground state at L = 
has been set to zero. 
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FIG. 3. Same as Fig. 2 for the exciton at v = 2/5 and 3/7. 
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FIG. 4. AV + T with N = 20 and T = 0.064e 2 /e/ (filled circles); and AV P (filled triangles) 
for 1/3 CF state. AV P (AV) is the interaction energy of the projected (unprojected) CF state. 
The statistical errors of Monte-Carlo estimation of the energies are shown on each point for the 
projected energy; in the case of the unprojected energies the error is of the order of the size of the 
filled circles. 
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FIG. 5. Estimation of AV P at the minimum of the exciton dispersion at v = 1/3. For each value 
of N, the minimum was determined by smooth polynomial interpolation in the neighborhood of 
the minimum to reduce the sensitivity of our thermodynamic estimation to the discreteness of L for 
finite systems. A linear x 2 -fit gives a thermodynamic estimate of 0.063 e 2 /elo, with an uncertainty 
of 2 (3) in the last digit if four (three) largest systems are used. 
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FIG. 6. The exact energies AV ex (dashes) and AV + T (filled circles) are shown for the exciton 
of the 1/3 state for systems of 6, 7 and 9 particles. T is chosen so as to match the two energies at 



the minimum. The exact energies in (c) are taken from Fano et al. 27]. 
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FIG. 7. Same as Fig. 4 for the 2/5 and 3/7 state. The exact energies in (a) are taken from 



d'Ambrumenil and Morf [29] and in (b) from He et al. [19] 
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FIG. 8. AV + T is shown for (a) v = 1/3, TV = 20, T = 0.064e 2 /d ; (b) 1/ = 2/5, iV = 30, 
r = 0.037e 2 /e/ ; and (c) i> = 3/7, iV = 27, T = 0.028e 2 /el . The typical error bar is shown on the 
first point in each case. 
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FIG. 9. for lowest energy exciton of (a) v = 1, N = 20, (b) v = 2, N = 30, and (c) v = 3, 
iV = 27. These correspond to the FQHE systems in Fig. 8. There are n minima/inflection points 
in the lowest branch of the collective mode of the v = n IQHE state. 
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FIG. 10. The charge density of the (a) v = 1 quasihole (solid line), which is the fully occupied 
lowest LL except for a hole at the origin, and (b) v = 1 quasielectron (dashed line), which is the 
fully occupied lowest LL plus an electron in the second LL at the origin. 
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FIG. 11. Exact A^T for (a) i/ = 1/3, N = 7; (b) v = 2/5, iV = 8. 
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FIG. 12. Monte Carlo estimate for AK using the unprojected CF exciton wave function for 
(a) 1/3 state of 20 particles, (b) 2/5 state of 30 particles, and (c) 3/7 state of 27 particles. The 
statistical error is smaller than the size of the symbol used. 
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FIG. 13. The kinetic energy of the 'unprotected SMA state', in units of huj c , for the collective 
mode of the 1/3 state (N = 20), shown by squares. The solid line plots Eq. (|2~8|) in unit of hu*. 



40 



